Critical wave-packet dynamics in the power-law bond disordered Anderson Model 
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We investigate the wave-packet dynamics of the power-law bond disordered one-dimensional An- 
derson model with hopping amplitudes decreasing as H nm oc \n — m\~ a . We consider the critical 
case (a = 1). Using an exact diagonalization scheme on finite chains, we compute the participa- 
tion moments of all stationary energy eigenstates as well as the spreading of an initially localized 
wave-packet. The eigenstates multifractality is characterized by the set of fractal dimensions of 
the participation moments. The wave-packet shows a diffusive-like spread developing a power-law 
tail and achieves a stationary non-uniform profile after reflecting at the chain boundaries. As a 
consequence, the time-dependent participation moments exhibit two distinct scaling regimes. We 
formulate a finite-size scaling hypothesis for the participation moments relating their scaling expo- 
nents to the ones governing the return probability and wave-function power-law decays. 

PACS numbers: 63.50.+X, 63.22.+m, 62.30.+d 



I. INTRODUCTION 

Non-interacting electron systems with uniformly dis- 
tributed disorder usually show an Anderson transition 
from localized to extended states. In general, such transi- 
tion occurs only for spatial dimensions greater than d = 2 
in the case of systems with short-range hopping, a result 
supported by a single parameter scaling theory QQ- On 
the other hand, when long-range couplings are assumed, 
a transition from localized to delocalized electronic states 
can be found even in ID disordered systems 0.0- In this 
case, one has an interplay between the hopping range and 
the degree of disorder. The former favors propagation 
while the later inhibits it. It is worthwhile to mention 
that propagation of carriers was also obtained in low- 
dimensional models with short-range hopping but pre- 
senting correlated disorder, such as random dimer chains 
ULEHjHH an d in chains with scale- free disorder ITol lll| . 
as well as in chains containing quasi-periodic structures, 
as for example Fibonacci, Thue-Morse, Harper sequences 

[3E0. 

On the other hand, for an ordered ID system with 
hopping terms decaying with a power-law with exponent 
a, it was shown that for a > 2 one recovers the result 
for the tight-binding model jl5| . More interesting is the 
behavior corresponding to < a < 1. For a — an 
initially localized wave-packet presents self trapping, i.e., 
the particle performs oscillations in a definite region of 
the lattice, visiting periodically the starting position. By 
increasing a the localization is lost. When the power 
exponent equals unity, and for sufficient short times, the 
packet diffuses with a diffusion coefficient that increases 
with the number of sites. This effect is absent in the 
model with only nearest-neighbor hopping [TBI . 

More recently, it was studied the dynamics of an elec- 
tron in a one-dimensional Anderson model with non- 
random hoppings falling off as some power a of the dis- 
tance between sites 16] . It was found that the larger the 



hopping range, the more extended the wave-packet be- 
comes as time evolves. When the disorder is increased, 
the wave-packet tends to remain more localized. For a 
low degree of disorder, the exponent a = 1.5 indicates 
the onset for fast propagation. Moreover, the inclusion 
of a dc electric field introduces the effect of dynamical 
localization. The fast propagation found for a < 1.5 is 
in agreement with the reported derealization of states 
located close to one of the band edges 0, 0, . 

The power-law random band matrix (PRBM) model 
also exhibits a delocalization transition 0,0. This model 
describes one-dimensional electronic systems with ran- 
dom long-range hopping amplitudes with standard devi- 
ation decaying as l/r a for sites at a distance r » b, 
where b is a typical bandwidth. It was shown that at 
a = 1 it presents an Anderson-like transition with all 
states being localized for a > 1 and extended for a < 1. 
At the critical point a = 1, the inverse participation ratio 
distribution, the wave- functions multifractal spectra, the 
level statistics and the time-evolution of the wave-packet 
size have been inv estig ated both analytically and numer- 
ically (lli lMEHL l2 a l23| l . Within the same spirit of the 
PRBM, a model for non-interacting electrons in a 2D lat- 
tice with random on-site potentials and random power- 
law decaying transfer terms was numerically investigated 
by exploring the finite-size scaling properties of the fluc- 
tuations in the mean level spacing 24] . It was found that 
the one-electron eigenstates become extended for transfer 
terms decaying slower than 1/r 2 . Finally, the Anderson 
transition in a ID chain with random power-law decay- 
ing hopping terms and non-random on-site energies was 
numerically characterized in ref. |25j. 

The moments of the position and those of the related 
probability density are known to exhibit different scaling 
behaviors [2(j. This feature reflects the complexity of 
the scaling laws governing the dynamics of quantum sys- 
tems [S]) . In this work, we investigate the critical dynam- 
ics of the power-law bond disordered Anderson model 
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where W nm is a random variable with a uniform dis- 
tribution in the interval [— 1,+1]. As a function of the 
exponent a characterizing the decay of the hopping am- 
plitudes, this model displays a localization-delocalization 
transition at a = 1 |25j in close connection with the 
PRBM [lEJ. For off-diagonal terms decaying slower than 
l/\n — m|, i.e., for a < 1, all states become delocalized. 
At a = 1 the states are critical and the level spacing 
statistics is between the Poison and Wigner surmises. In 
the limit of a — > oo one recovers the ID Anderson model 
with just first-neighbors random hopping amplitudes. 



FIG. 1: The average participation moment (P q ) as a function 
of the chain size (TV) for several values of q ranging from 2 to 
6. The distinct fractal dimensions (P q oc N "^'^) exhibited 
by the participation moments reflect the multifractality of the 
critical eigenstates. 



(a = 1) by looking at the moments of the wave-function 
probability distribution. Using an exact diagonalization 
scheme on finite chains, we compute the participation 
moments of all energy eigenstates and follow the time- 
evolution of an initially localized wave packet. After that, 
we perform an average of the participation moments over 
different configurations and energies, which allow us to 
compute the critical exponents D q associated with the 
multifractal character of the critical eigenfunctions. Ex- 
amining the time-evolution of an initially localized wave- 
packet, we obtain the decay exponent of the autocorre- 
lation function (C(t) oc t~ D2 ) and the size dependence 
of the asymptotic return probability (R(t) oc N~ D z ). 
We will employ a finite-size scaling analysis of the time 
dependent participation number moments, relating their 
scaling exponents with the power-law exponents of the 
evolving wave-packet. 



II. MODEL, FORMALISM AND RESULTS 

A. Model Hamiltonian 

We consider a single electron in a ID chain with open 
boundaries, described by the Anderson Hamiltonian 
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H 



|n) (m\ 



(i) 



where |n) represents the state with the electron localized 
at site n. In the present pure random bond Anderson 
model, the on-site potentials are site independent and 
in Eq. 1 were taken to be ej = without any loss of gener- 
ality. Long-range disorder is introduced by assuming the 
hopping amplitudes Hij to be randomly distributed and 
also displaying a power-law decay. We will, hereafter, 
consider 



(2) 



B. Eigen- functions participation moments 

We will be particularly interested in investigating the 
participation moments and the time-evolution of an ini- 
tially localized wave-packet at the critical point a = 1. 
The participation moments for a particular disorder con- 
figuration and eigenstate will be defined as the inverse of 
the moments of the probability density 
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where fn is the amplitude at site n of the j-th eigenstate 
obtained from an exact diagonalization scheme on finite 
chains with sizes ranging from N = 400 up to N = 3200 
sites. Here we used 32 x 10 3 states for each chain size 
in order to average over distinct disorder configurations 
and energies 
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where <> stands for configuration average. Therefore, 
P q probes the average q-th. moment of the critical sta- 
tionary states. In fig. 1 we show the finite-size scaling 
of different moments ranging from q = 2 up to q = 6. 
As it has been shown by Wegner 28], all moments show 
a power-law dependence on the form (P q oc N Dqt - q ~ 1 ' > ). 
Each moment has its characteristic exponent D q as indi- 
cated in Fig. (1) which reflects the multifractality of the 
critical eigenstates. Our results indicate that the fractal 
dimension D q slowly decreases from D2 = 0.676 towards 
the asymptotic value Drx, = 0.5. 



C. Wave-packet Dynamics 

In order to obtain the time-evolution of an initially 
localized wave-packet (|$(£ = 0))), we expand the wave- 
function in the Wannier representation 
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FIG. 2: (a)The return probability function for t — > oo as 
a function of the number of sites. The small power-law de- 
cay exponent R(t -> oo) oc N~ D 2 , with Z)f = 0.607(10), is 
consistent with a wave-function with a slowly decaying enve- 
lope, (b) Autocorrelation function versus time t for chains 
with N — 400 up to N = 3200 sites. The power-law decay for 
intermediate times with an exponent D2 = 0.676 is in good 
agreement with the previous numerical estimate for Z?2. 



We solve the time-dependent Schrodinger equation for 
the wave-function components f n {t) {% = 1) 



df n (t) 

dt 



N 



:/»(<). 



(6) 



using the numerical formalism proposed in ref. |29| : 

= U f exp (-iDt)U|$(t = 0)), (7) 

where D is the diagonal form of the Hamiltonian and 
U is a unitary matrix. In what follows we consider the 
electron initially localized in a single site, i.e. = 
0)) = I no). We start by discussing the behavior of the 
autocorrelation function: 



C(t) = - f R(t)dt, 
t Jo 



(8) 



where R(t) = \f no (t)\ 2 denotes the return probability. 
The time dependent autocorrelation function provides in- 
formation about localization of the wave-packet as well 




FIG. 3: Average wave-function profile after a very long evo- 
lution time (t — > 00). In the inset, we note that the envelope 
displays a power-law decay with an exponent 1 — D* = 0.393, 
is agreement with the measure exponent of the asymptotic re- 
turn probability. 



as of its fractal character. In fact, after an initial waiting 
time, the autocorrelation function shall behave as: 



C(t) oc r 



(9) 



for one-dimensional systems, which be definition is the 
same power-law decay presented by the return probabil- 
ity. In fig. 2(a) we plot the return probability for t — > 00 
as a function of the number of sites. After the initial 
decay due to the wave-function spread, it saturates at 
a size-dependent plateau which decreases as 1/N D2 . In 
fig. 2(b), the autocorrelation function versus time t is 
shown for different chain sizes. It exhibits a power depen- 
dence for intermediate times whose exponent is in good 
agreement with the previous estimate for D%. Its satura- 
tion for t — > 00 has the same finite-size origin presented 
by the return probability. 

The asymptotic return probability decaying slower 
than 1 /N indicates that the statistically stationary wave- 
function does not assume a uniform profile. In fig. (3) 
we plot the average stationary wave-function profile 
achieved long after the wave-function reflection at the 
chain boundaries. It develops a power-law decaying tail 
with \f n (t — ► oo)| 2 oc l/(n — rio) 1_£>2 (see inset), as 
expected according to the wave-function normalization 
constraint. Before reaching the chain boundaries, the 
average wave-function envelope can be well described by 



\f(n,t)\ 2 =R(t)(n-n + l)-^- D ^ 



(10) 



for n — no < N*(t), where N*(t) is a time depen- 
dent cutoff after which the wave-function decays expo- 
nentially and, therefore, can be considered as having a 
vanishing amplitude. The return probability decays as 
R(t) oc t~ D2 . It is straightforward to show, using the 
normalization criterium, that the cutoff shall scale as 
N*(t) oc t D2 /- D 2 , The above asymptotic power-law shape 
of the wave-function is in agreement with the previous 
prediction for the wave-packet spreading in quantum sys- 
tems with fractal energy spectra and eigenfunctions|27j| . 
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FIG. 4: The normalized average participation number P2/N 

as a function of the scaled time (t/N D2 ^ D ' 2 ). After a short 
transient time, the data collapse indicates that the participa- 
tion number grows as t £>2 / D 2 anc j saturates after the wave- 
function reflection at the chain boundaries. The destructive 
interference which sets up during the reflection at the bound- 
aries is signaled by a maximum in P 2 . The data collapse 
corroborates the scaling hypothesis in Eq. (12). 
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We further calculated the time dependent participation 
moments 



PS) = 



1 
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(11) 



The second moment (P2) has been commonly used as a 
measurement of the number of sites that participate in 
the wave function. For a wave-packet fully localized at 
a single site, P 2 — 1, while for a uniform distribution of 
the wave function, P2 reaches the maximum value TV. In 
general, P2 is size independent in the localized phase and 
scales linearly with the system size for extended states. 

In the present case, for which the wave-function de- 
velops the power-law form given by Eq.(10), the par- 
ticipation moments present a non-trivial scaling. Two 
distinct behaviors can be anticipated. For q < 1/(1 — 
D2), the participation moment shall display a power- 
law growth (after a short waiting time) on the form 
P q (t) oc t^ -1 )- 02 / 1 ^ i followed by a saturation after reach- 
ing the chain boundaries. The saturation plateau is pre- 
dicted to scale as P{N,t -> 00) oc N^l Both time 
regimes can be represented by a single finite-size scaling 
hypothesis for the participation moment which can be 
written as 



P q (t,N) = N^ q -^g(t/N D ^ D2 ) 



(12) 



with g(oo) = cte and g(x -> 0) oc x (<2-i)£>2/.d* _ In fig 4j 
we show our numerical results for the scaled second mo- 
ment (P2/N) versus scaled time (t/N D ^ /° 2 ), which falls 
in the above regime of q < 1/(1 - D$). These data were; 
averaged over 100 time histories considering distinct dis- 
order configurations. The second moment remains con- 
stant during a very short initial waiting time before the 
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FIG. 5: Scaled participation moments P q (N,t)/N qD z as a 



function of the scaled time (t/N 2 



for q = 3, 4, 5. The 



data collapse corroborates the scaling hypothesis in Eq. (15). 
The maximum, also appearing in fig. 4, signals the destructive 
interference between the incident and reflected waves near the 
chain boundaries. 



spread effectively starts to take place. After that, we 
found that the above scaling hypothesis reproduces very 
well the participation time-evolution as data from dis- 
tinct chain sizes could be collapsed in the same curve. 
As anticipated, the participation moment has a power- 
law growth followed by a saturation. The slight devia- 
tion from perfect collapse is due to corrections to scaling 
which are more pronounced for the smaller system sizes. 
The maximum displayed by the participation number re- 
flects the destructive interference between the incident 
and reflected waves near the chain boundaries. 
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For q > 1/(1 — £>*), the presence of a cutoff does not 
influence the initial spread of the participation moments. 
Therefore, the initial growth is completely determined by 
the return probability decay and results in P q (t) oc t qD2 . 
On the other hand, the scaling of the asymptotic plateau 
only depends on the wave- function decay exponent, being 
given by P(N, t — > oo) oc N qD ^ . For this regime, the 
proper scaling form of the participation moment reads 

P q (t,N) = N qD * g'(t/N D ^ D2 ) (13) 

with g'(oo) = cte and g'{x — > 0) oc x qL>2 . In figure 5 we 
show the scaling analysis of participation moments for 
this regime {q — 3,4, and 5). The trends are mainly the 
same as that presented by Pi , namely data collapse with 
saturation for long times and a power-law growth for in- 
termediate times preceded by a size independent initial 
waiting time. However, the participation moment scal- 
ing exponent now depends on the wave-function decay 
exponent as predicted by Eq.(13). 

In both regimes, the characteristic time scale after 
that finite-size effects starts to take place scales as t* oc 
n d%/d 3) with Df/D 2 = 0.90(2). The fact that this scal- 
ing exponent is smaller than unity (and therefore corre- 
sponds to a faster than ballistic dynamics) is related to 
the non-local nature of the long-range couplings. Finally, 
for exponentially localized wave- functions [-Pg] 1/ ' ^ is 
proportional to the localization length for any q. Using 
N* as a measure of the characteristic size of the present 
power-law decaying wave-packet, this feature is only true 
for q < 1/(1 — D%), For higher moments it displays a 
sub-linear size dependence whose exponent continuously 
decreases towards £)* as q — > oo. 

III. CONCLUSIONS 

We investigated the one-dimensional Anderson model 
with off-diagonal disorder and matrix elements Hij de- 
caying as l/\i — j\ a for a — 1 (critical regime). Using an 
exact diagonalization scheme on finite chains, we com- 
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